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Abstract 

Degeneracies in the spectrum of an adiabatically transported quantum system are important to 
determine the geometrical phase factor, and may be interpreted as magnetic monopoles. We inves- 
tigate the mechanism by which constraints acting on the system, related to local symmetries, can 
create arbitrarily large monopole charges. These charges arc associated with different geometries of 
the degeneracy. An explicit method to compute the charge as well as several illustrating examples 
are given. 
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1 Introduction 



When a quantum system is adiabatically transported along a closed loop C in a parameter space, in 
addition to the usual dynamical phase every non-degenerate eigenstate accumulates an extra phase A<1>. 
For a 3D-parameter space, M. V. Berry gives a geometrical interpretation of A<3? as the flux of a 
magnetic- like field B through a surface with boundary C (Berry, 1984). The degeneracies of the 
spectrum in parameter space are the singularities of the vector field B, and therefore play an important 
role in connexion with the geometric phase. Each degeneracy can be seen as a charge distribution 
located at the contact point between energy surfaces. Because the eigenstates are smooth and single 
valued outside the degeneracies, the total charge of the distribution, i.e. the monopole charge or 
moment, is necessarily an integer multiple of the elementary charge go = 1/2. In the generic case of 
a diabolical contact (Berry, 1984), the monopole charges are precisely ±<?o- However, higher integer 
multiples of go may occur. For instance, for light propagating through a twisted anisotropic dielectric 
medium there are experimental situations (Berry, 1986) where the monopole charges are ±2go- Other 
examples arise in condensed matter physics. In a bidimensional periodic crystalline lattice, the Hall 
conductance of a gas of independent electrons is proportional to a topological index, the Chern index, 
that measures the net charge inside a closed surface associated to the first Brillouin zone (Simon, 1983). 
In some models, large jumps of the Chern index that cannot be explained by elementary charges have 
been observed (Leboeuf et al., 1990) (Faure & Leboeuf, 1993). 

Our purpose is to discuss a generic mechanism for the production of monopole charges larger than 
go- The mechanism is due to constraints that act on the system, and may be associated to local 
symmetries. Section El introduces the notation and some general formulae for A$. The generic case 
of a diabolical contact is briefly recalled in section |3 In section 01 which is the central part of the 
paper, a broader class of situations that incorporates constraints is considered. The consequences 
on the geometry of the contact point and on the corresponding monopole charges are discussed, and 
a method to explicitly compute them is introduced. The latter provides an alternative and simple 
view of the quantization of the monopole charge as a sum over winding numbers associated to Dirac 
strings. We illustrate the general approach by two examples in section El an d show how arbitrarily 
large multiples of go can actually be generated in [HJ The last section contains a discussion of some of 
the results obtained. 

2 Background material 

We consider a quantum system governed by a Hamiltonian depending on three real parameters r = 
(x,y, z). We denote V the parameter space and suppose that in the neighborhood of r = two 
eigenvalues are close to each other. Hence, the coupling to other states can be neglected and the 
Hamiltonian can be restricted to a bidimensional eigenspace. The most general form of the reduced 
Hamiltonian is 



s(r) and the three components of e(r) = [e x (r),e y (r),e z (r)] are smooth real functions in "P. er = 
[a x ,a y ,a z ] are the Pauli matrices. 

The eigenvalues of H (r) are given by 
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(2) 
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where e = w + + e|. The points r of the parameter space where degeneracies occur are given by 

the set M where e vanishes. 

One possible orthonormal eigenbasis corresponding to E±(r) is 

l^(r)> d = f ; ( ex(r) ~ iey(r) ) ■ (3) 

2e(r) (e(r)=Fe*(r)) \ ±e ^) ~ e z (r) J 

It follows from this expression that the eigenstates of H are not only singular on Ai but more generally 
on a larger subset T> of V . With the choice (jHJ), T> is given by the points where e(r) =p e z (r) vanishes. 
For convenience, we split P into two sets, T> + and X>_, that corresponds to e x (r) = e y (r) = and 
e z (r) > (e z (r) < 0), respectively. These two sets intersect on M.. 

Let C be any closed loop in V not intersecting V. Berry (Berry, 1984) gave a geometrical in- 
terpretation (that is, coordinate free in V) of the non-dynamical part A<I>-|- of the phase shift that 
\i[>±(y)) acquires when H(r) follows C adiabatically. Simon (Simon, 1983) completed the picture with 
a topological interpretation of Berry's phase in terms of connection on a suitable fiber bundle. Topo- 
logical arguments were anticipated in (Stone, 1976, see in particular the section: "A topological test 
for intersections", p. 85). The phase shift after one traversal of C is given by the circulation 

A$+ = -A$_ = A3> d = / A(r). all (4) 



c 



of the vector A = f — Im 



(^+|v r |v>+> 



. By Stokes's theorem, A<3? can alternatively be computed from, 
A$ = / B(r).dE . (5) 
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S is any surface in V not intersecting T> whose border is C. B(r) is the 3- vector field (more generally a 
2-form) obtained by taking the curl of A. Specifically, outside T> and in the two-dimensional subspace 
where our analysis is restricted 

B = {e + -e_y • (6) 



From (J2), (j2J and (J3J) we may write, outside V, 



A = e v Vr ^ e ^ re y (i + 2i \ ) (7) 



and ^ 
B= 2^ 



e x (V r e y A V r e z ) + A V r e x ) + e 2 (V r e x A V r e y ) . (8) 



3 Generic case (no constraints) 

Generically M is either empty or an isolated point (see (von Neumann & Wigner, 1929)). We consider 
the latter case and work in a local chart centered at the contact point, that we assume is located 
at r = 0. Since generically 

A = f det[V r e x , V r e y , V r e 2 ]| r=Q (9) 
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is different from zero, by virtue of the local inversion theorem we can consider (e x , e y ,e z ) to be a new 
local chart in a neighborhood of r = 0. In other words, we can take r = e. This is the generic case. 
The two energy surfaces E±(r) intersect through an isolated doubly-conical contact point (diabolical 
point) (von Neumann & Wigner, 1929). From (JSJ) we have 



Close to a conical intersection the vector B corresponds therefore to the field created by a magnetic 
monopole of charge g = go = 1/2 (Dirac, 1931; Dirac, 1978) (for the state \ip-(r)) the opposite charge 
must be taken). A(r) is the vector potential associated to it. The sets T>± are the well known Dirac 
(half) strings which physically can be interpreted as semi-infinite solenoids (with end points at r = 0) 
carrying a magnetic flux ±47r<70) respectively. This magnetic analogy cannot however be extended 
globally since there is no superposition theorem available here. 

4 Finite number of constraints 
4.1 Non-elementary monopoles 

We now consider situations where some constraints acting on the system are present. These con- 
straints can be thought of as independent real control parameters A = {^n}ne{o,...,C-i} that impose C 
independent relations between the derivatives of e at r = 0. Physically the control parameters are 
distinct from the adiabatic parameters r in that they reflect a local symmetry in V and therefore do 
not change during the adiabatic transport along C. Due to these additional constraints, it may hap- 
pen that the determinant A defined by Eq. © vanishes. The latter condition requires generically the 
existence of one external parameter A in addition to r. In that case, the contact between the energy 
surfaces restricted to the r -space V is not a diabolical point anymore (see figure For instance, 
the charge g = 2go occasionally observed in Ref. (Leboeuf et al., 1990) was found to correspond to a 
parabolic contact. By increasing the number of parameters A one can accordingly cancel an increasing 
number of derivatives of e at r = 0, and therefore increase accordingly the degree of tangency of the 
surfaces defined by ©. 

If the number of constraints C is finite, in the full (C + 3)-dimensional parameter space (that 
includes both the adiabatic and the control parameters), the set of degeneracies is C-dimensional and 
generically crosses transversely at an isolated point located at r = any 3-dimensional space defined 
by A = constant. It requires an infinite number of constraints to transform A4 into a continuous family 
of points since an infinite number of vanishing derivatives of e is needed. This is the situation when, 
for instance, the system is invariant under a global symmetry, like the time reversal symmetry. In that 
case e y = at any point in V and M is a 1-dimensional sub-manifold of V . Global symmetries that 
increase the dimension of the contact manifold are not considered here, we restrict to the cases where 
the number of constraints due to local symmetries is finite. 

Under these assumptions, it follows from (jHJ that divB is a distribution with support in r = 0. It 
can therefore be expanded as, 



As for an arbitrary circuit in parameter space (Berry, 1984, p. 55), and in contrast to the diabolical 
contact, generically in the two-dimensional subspace of the degeneracy the geometric phase is not 
completely determined by (|11|) because of the non- vanishing curl of (jHJ). Our purpose is to determine 




(10) 



divB = gS(r) + terms involving derivatives of 5(r). 



(11) 
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order of tangency 

Figure 1: The unfolding shown in this picture illustrates the changes in the geometry of the contact 
point between two energy surfaces induced by the variation of some control parameter A. 



an explicit rule to compute g. A more detailed qualitative discussion of the contribution to the 
geometric phase of the curl of B and of the additional multipolar terms in (jllj) will be made in the 
final section. 

4.2 The monopole charge as a sum of winding numbers 

One of the most remarkable properties of the monopoles is the quantization of their charge, which can 
only take integer multiple values of go = 1/2. This easily follows from a topological argument given 
by (Stone, 1976). Consider a small loop C that is smoothly retracted to a point without crossing D. 
Because the eigenstates remain smooth and single valued, A$ must tend to an integer multiple of 2n. 
But A<1? is the flux of B through any surface S with boundary C. When C shrinks to a point, 
S tends either to a point or to a closed finite surface that possibly encloses a monopole charge g. 
Therefore A<E> — > 4-7r<7, which must be an integer multiple of 2tt. The quantization of g follows. 

This argument is, however, unable to provide a method to compute the charge. One usually has 
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to explicitly compute the integral © (using (jSJ for the magnetic field) over a sphere enclosing the 
degeneracy. In the case of a conical contact, the structure of the integral is transparent: B. cffi is 
half the solid angle subtended from the origin by cffi. However, for an arbitrary contact point in 
the presence of additional constraints B.dS is a complicated function of the polar and azimuthal 
angles. The structure of the two-dimensional integral is therefore non-trivial and in general difficult 
to compute. It obscures moreover the simplicity of the result (an integer multiple of 2ir). 

Rather than the magnetic field B, from now on we work with the potential A. As we will show, the 
2D-flux integral across the sphere is replaced by one or several ID-circulation integrals. This scheme 
leads to a general and simple rule for calculating the charge that explicitly exhibits the quantization. 

When the number of constraints C is finite, the subset in V defined by e x (r) = e y (r) = is 
the union of the algebraic curves defined near r = by the first non-vanishing terms of the Taylor 
expansion of e. The expected "pathologies" of this subset are intersections of the curves at r = 0. 
When the conditions e z (r) > and e z (r) < - that define the Dirac strings T> + and T>- - are added, 
the number of curves can possibly be reduced by half. In a small enough neighborhood of r = 0, 
T> + and 2?_ are made of sets of half curves starting at the origin. Locally one cannot determine if a 
given couple of half strings belongs to the same algebraic curve or not. We also stress that it is the 
choice of gauge which fixes T> + . Adding a gradient to A will not change the value of A$ but can 
modify T> + significantly by changing the position and the number of half strings. What really matters 
is the algebraic sum of the flux they carry. To illustrate this point, a useful gauge transformation is to 
make a rotation R £ SO(3) on e: If U denotes a unitary representation of SO(3) in the bidimensional 
Hilbert subspace, the relation 



together with implies that such a rotation corresponds to a unitary transform of the eigenvec- 
tors of H and hence will not affect the phases. Yet a simple permutation of the components of e 
can change T> + drastically. For example, a rotation by ir about e x changes the sign of e z without 
modifying e x and e^. One can therefore exchange T> + and It is moreover not excluded that a 
gauge transformation that considerably simplifies the sets T) + and T>- exists in general, by reducing 
them for example to a single half Dirac string. But no general rule for constructing such a gauge 
transformation is available. 

From now on we fix the gauge according to Eq. (JJJ). S is chosen to be a surface which is diffeomor- 
phic to a sphere centered at r = but pierced with N holes at the level of the strings contained in T> + . 
We denote Ci the (oriented) boundary of the Ith hole. It is a small loop of typical radius 5i > that 
encircles the ith half string. Then the integration contour in @ reduces simply to dS = C = 5^/3= i A 
(see figure El • 

Since C can be taken arbitrarily close to T> + , and not intersecting we have 




(12) 



e(r) 



1 + 0(Si) when re£ ( . 



(13) 



From this it follows that 




lim — — 

closed sphere 47T 
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(14) 



i=i 



where 




(15) 
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Figure 2: Example of a set T> + in the presence of constraints. On the left part of the figure is displayed 
a closed Dirac string crossing itself once at r = 0, in addition to a half Dirac string. From a global 
perspective, one can predict that the former will not contribute to the monopole charge. Nevertheless, 
in a neighborhood of r = (right figure) one has to compute the winding of Z\ when r traverses Ci 
for each of the five half-strings. 



wi has a natural topological interpretation. It is the algebraic winding number of the complex number 

Z z (r) = e:r (r)+ie y (r) (16) 

when r traverses C\. It is a well defined quantity: Zj(r) never crosses zero because r never touches T>. 

Eq. (|14jl yields an alternative interpretation of the charge as a sum of winding numbers associated 
to the Dirac strings. It also provides a method to compute the charge, or to generate arbitrary values. 

5 Two examples 

We end up by illustrating the procedure with several examples (the generic case of section|3]is similar to 
the second one). In all the cases described below, we will Taylor expand e near the origin and suppose 
that the terms that we explicitly retain are sufficient to fix the local geometry of T> completely. 
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5.1 An example with two choices of gauge 

Consider cl CELS6 where in a suitable choice of coordinates we have, 

e(x,y,z) = [x 2 + z, y 2 + z, z] + higher order terms . (17) 

It follows immediately that T> + = {0} so that g = 0. To illustrate the importance of the gauge, we now 
treat the same example but in a gauge that exchanges T) + and 2?_ . T>- is made of four half branches 
of parabolae: {[x = z i V = ±\/ — z ) z — 0]}- Close enough to r = they can be assimilated to 
their tangent, i.e. the two axes {x = ±y,z = 0}. The surface S looks topologically like the one in 
figure El 




Figure 3: Near r = 0, the Dirac strings in 2?_ are tangent to the two straight lines x = ±y, and the 
four loops Ci are circles of radius S and constant y with center in one of the lines. 

The four small loops Ci can be chosen to be centered on P_ and parametrized by [x = ±e + 
5 cos 8, y = ±e, z = —e 2 — sgn(y)5sin#]. 9 runs from to 2ir, whereas e and 5 are fixed and small 
strictly positive quantities measuring the distance to the origin and the radius of the loops, respectively. 
From the definitions ((TB|) and (fTT|) we get Z = 5[— sgn(y)(l + i)sin# + 2 sgn(x)ecos 6] + 0(5 2 ). Z 
describes (a small perturbation of) an ellipse centered at the origin whose great axis is the segment 
[—(5(1 + i), (5(1 + i)] and whose width is proportional to e. The corresponding winding number is 
sgn(x)sgn(y) and therefore the total charge given by (|14|) vanishes. 
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5.2 An example with a finite but arbitrarily large number of constraints 

Now consider the case 

e(x,y,z) = [(xy) n , ^(x 2n — y 2n ), z] + higher order terms , (18) 

where n is a strictly positive integer. The propagation of light in a twisted anisotropic dielectric 
medium considered in (Berry, 1986, see eq. (45)) corresponds to n = 1. T> + is the half-axis x = y = 
0, z > 0. There is only one loop to consider (see figure @J). 




Figure 4: The examples considered in section 15.21 and El lead to one half Dirac string x = y = 
and z > 0. 



Let us fix z = zq > and S > 0. Consider the loop 
C = [x = 5 cos 9,y = 5 sin 8, z 
Then A<I>/(27r) is the winding number around zero of 

'sin(2#)\ n 



9 £ [0,2vr[ 



Z{0)=8 



2n 



+ 2 ( cos 



2n 



sin 



2d 



(19) 



(20) 



If n is even there is no winding, since the real part of Z remains always positive. Therefore g = 0. 
In contrast, if n is odd 7.(6) winds clockwise twice, and therefore g = —2go. Figure |S] illustrates two 
different contact points corresponding to the two possible charges. 
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6 Arbitrarily large charges 



In order to get arbitrarily large winding numbers for Z, we make use of Chebyshev's polynomials, 
defined by (Magnus et al., 1966, chap.V, §5.7) 



Qi(<7) = cos(narccosi 



1 



(21) 



where n is an arbitrary integer and q a real number satisfying — 1 < q < 1. Let us take 



/ 



e(x,y,z) 



y/x 2 + y 2 



i(Vx 2 + y 2 ) C r 



COS 



/ 7T N 



X 



+ sin 



\2nJ ^ x 2 + yi \2nJ ^ x 2 + y 2 



+ higher order terms . 



V * / 

(22) 

f(q) and <?((/) are some regularizing real functions that vanish at q = faster than any singular 
denominator coming from C n . An appropriate choice is, for instance, f(q) = g(q) = q n . Though C n 
vanishes n times in [—1,1], one can always choose a small neighborhood where r = is the unique 
point in M. As in the last example, the set P+ is made of one Dirac half string on the half-axis x = 
y = 0, z > 0. Constructing the same loop as before one gets 

Z(0) = f{5) cos(nfl) + i g{5) sm(n6) , (23) 

which corresponds to g = ng$. This provides an explicit example of a monopole carrying an arbitrar- 
ily large charge, and therefore creating arbitrarily large phase shifts for a loop located close to the 
degeneracy in parameter space. The contact between energy surfaces with n = 2 is illustrated in part 
(a) of figure El which is identical to the geometry of the contact defined by Eq.(|18|) with n = 1 (up to 
an irrelevant scaling factor). 



7 Discussion 

Spectral degeneracies with charges higher than one have been observed in several contexts. For 
instance, via the variations of Chern numbers in quantum systems with mixed and chaotic classical 
dynamics (Leboeuf et al., 1990; Leboeuf et al., 1992). These integer topological invariants measure 
the transport (Hall conductivity) of each band of a doubly periodic system. Values of g/go = 2 and 3 
were found (cf in Fig. 2 of Ref. (Leboeuf et al., 1990) the contact points at 7 = 0.7106 between the 
second and third levels, and at 7 = 0.4477 between the fourth and fifth levels, respectively). A closer 
analysis of the geometry of the contact points reveals that these two degeneracies coincide with the 
parabolic and cubic contacts described by Ea. (|22j) with n = 2 and 3, respectively. For instance, the 
g/go = 2 contact is parabolic in the Bloch angles and linear in the parameter 7. 

We have also exhibited the possibility of having a contact point with total charge equal to zero. 
Although its presence will in general be difficult to detect, it would be interesting to find an explicit 
quantum mechanical model displaying this curious feature. 
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Figure 5: Geometry of the contact point between energy surfaces for z = 0. Part (a) corresponds 
to Ea. Q18|) with n = 1 and to Ea.([22|) with n = 2. The charge is g/go = 2. Part (b) corresponds 
to Eq.(|18jl with n = 2. The charge vanishes. Notice the difference in the geometry of the contact 
between both cases. 



In the electromagnetic analogy based on the general formula, the "sources" of the geometric phase 
in parameter space are the spectral degeneracies plus some additional currents not related to them. 
In the particular (and generic) case of a diabolical contact, the contribution of the degeneracy is 
simple: it acts as a pure monopole charge, thus contributing to divB but not to curlB. However, 
for an arbitrary contact point - defined by Eq.(|%|) - the situation is different. On the one hand, the 
divergence of B is not simply determined by a charge, but more generally by a charge distribution 
(cf Eq,(|ll[)). The latter can contain higher multipole moments. Moreover, the curl of B need not be 
zero. In the present study, we have only considered a particular aspect of the sources of geometric 
phase associated to the field produced by a degeneracy, namely the total charge of the distribution. 
A more complete classification scheme is clearly needed. 
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